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Abstract. We use combinatorial and Fourier analytic arguments 
to prove various non-existence results on systems of real and com- 
plex unbiased Hadamard matrices. In particular, we prove that 
a complete system of complex mutually unbiased Hadamard ma- 
trices (MUHs) in any dimension d cannot contain more than one 
real Hadamard matrix. We also give new proofs of several known 
structural results in low dimensions, for d < 6. 



1. Introduction 

A new approach to the problem of mutually unbiased bases (MUBs) 
was recently given in [IE], based on a general scheme in additive com- 
binatorics. In this paper we continue the investigations along this line, 
and prove several non-existence results concerning complete systems of 
MUBs, as well as some structural results in low dimensions. Let us 
remark here that the existence of MUBs is equivalent to the existence 
of mutually unbiased Hadamard matrices (MUHs) as explained below. 
In most of the paper it will be more convenient to deal with MUHs. 

The paper is organized as follows. The Introduction contains a 
standard summary of relevant notions and results concerning MUBs 
and MUHs. We also recall some elements of the general combinatorial 
scheme which was used in [16J. In Section [2] we use discrete Fourier 
analysis to prove several structural results on MUHs in low dimensions. 
Finally, in Section [3] we prove non-existence results including the main 
result of the paper: a complete system of MUHs can contain at most 
one real Hadamard matrix. We also give a new proof, without using 
computer algebra, of the fact the Fourier matrix F 6 cannot be part of 
a complete system of MUHs in dimension 6. 
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Recall that two orthonormal bases in C d , A = {ei, . . . , e^} and B = 

{fi, . . . , id} are called unbiased if for every 1 < j, k < d, \ (e,, i k ) \ = —j=. 

In general, we will say that two unit vectors u and v are unbiased if 

|(u, v)| = —=. A collection Bo, . . . B m of orthonormal bases is said to be 
yd 

(pairwise) mutually unbiased if every two of them are unbiased. What 
is the maximal number of pairwise mutually unbiased bases (MUBs) in 
C d ? This question originates from quantum information theory and has 
been investigated thoroughly over the past decades (see [9] for a recent 
comprehensive survey on MUBs). The following result is well-known 
(see e.g. [DEI [20]): ' 

Theorem 1.1. The maximal number of mutually unbiased bases in C d 
is at most d + 1. 

Another important result concerns prime-power dimensions (see e.g. 

mini Era). 

Theorem 1.2. A collection of d + 1 mutually unbiased bases (called 
a complete set of MUBs) exists if the dimension d is a prime or a 
prime-power. 

However, if the dimension d = p" 1 . . . p^ k is composite then very little 
is known except for the fact that there are at least pj J +1 mutually 
unbiased bases in C d where p" J is the smallest of the prime-power 
divisors. In some specific square dimensions there is a construction 
based on orthogonal Latin squares which yields more MUBs thanpj' + l 
(see [19J). It is also known [TS] that the maximal number of MUBs 
cannot be exactly d (i.e. it is either d + 1 or strictly less than d). 

The following basic problem remains open for all non-primepower 
dimensions: 

Problem 1.3. Does a complete set ofd+l mutually unbiased bases 
exist in C d if d is not a prime-power? 

The answer is not known even for d — 6, despite considerable efforts 
over the past few years ( [21 EJ [7J [13] ) . The case d = 6 is particularly 
tempting because it seems to be the simplest to handle with algebraic 
and numerical methods. As of now, numerical evidence suggests that 
the maximal number of MUBs for d = 6 is 3 (see [6j [7J El |2T]). 

It will also be important for us to recall that mutually unbiased 
bases are naturally related to mutually unbiased complex Hadamard 
matrices. Indeed, if the bases B , . . . , B m are mutually unbiased we 
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may identify each B\ = {ef , . . . , e^} with the unitary matrix 

rrri [/„(o) 
i u i\j,k — \ ■? ' k 

i. e. the k-th column of Ui consists of the coordinates of the k-th vector 
of Bi in the basis Bo. (Throughout the paper the scalar product (., .) 
of is conjugate-linear in the first variable and linear in the second.) 
With this convention, Uq = I the identity matrix, and all other matri- 
ces are unitary and have all entries of modulus l/y/d. Therefore, for 
1 < I < m the matrices Hi = \fd~Ui have all entries of modulus 1 and 
complex orthogonal rows (and columns). Such matrices are called com- 
plex Hadamard matrices. It is thus clear that the existence of a family of 
m + 1 mutually unbiased bases Bo, ... , B m is equivalent to the existence 
of a family of m complex Hadamard matrices Hi , . . . , H m such that for 
all 1 < j ' k < m, -^H*Hk is again a complex Hadamard matrix. 
In such a case we will say that these complex Hadamard matrices are 
mutually unbiased (MUHs). A system Hi, ... , H m of MUHs is called 
complete if m = d (cf. Theorem 11.11) . We remark that there has been 
a recent interest in real unbiased Hadamard matrices jH [Til EE] > an d 
the main result of this paper is that no pair of real unbiased Hadamard 
matrices can be part of a complete system of MUHs (see Corollary 
3.21) . The system Hi, . . . H m of MUHs will be called normalized if the 
first column of Hi has all coordinates 1, and all the columns in all the 
matrices have first coordinate 1. It is clear that this can be achieved 
by appropriate multiplication of the rows and columns by umimodular 
complex numbers. We will also use the standard definition that two 
complex Hadamard matrices Hi and H 2 are equivalent, Hi = H 2 , if 
Hi = D1P1H2P2D2 with unitary diagonal matrices Di, D 2 and permu- 
tation matrices Pi,P 2 - 



The crucial observation in [16] is that the columns of Hi, . . . ,H m 
can be regarded as elements of the group Q = T d , where T stands 
for the complex unit circle. By doing so, we can use Fourier analysis 
on Q to investigate the problem of MUHs. We will now collect some 
notations that will be used in later sections (the notations in this paper 
are somewhat different and more convenient than in [IS]). The group 
operation in Q is complex multiplication in each coordinate. The dual 
group is Q = Z d , and the action of a character 7 = (n, r 2 , . . . , r d ) 6 7h d 
on a group element v = {vi,v 2 , . . . ,Vd) € T d is given by exponentiation 
in each coordinate 7(v) = v 7 = v^v^ 2 . . .v r d d . The Fourier transform 
of (the indicator function of) a set S C Q is given as S(j) = J2 s es s7 - 
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As in [16], introduce the orthogonality set ORTd = {v = (vx, . . . , v d ) £ 
T d : t> i + ■ • • + v d = 0}, and the unbiasedness set UB d = {v = 
(vi, . . . ,Vd) G T d : |t>! + • • • + t> d | 2 — d = 0}. Then the (coordinate- 
wise) quotient v/u of any two columns from the matrices Hx, ■ ■ ■ H m 
will fall into either ORTd (if v and u are in the same matrix) or into 
UBd (if v and u are in different matrices). This enables one to invoke 
the general combinatorial scheme which we called "Delsarte's method": 
we refer the reader to (THj for the details. 

2. Structural results on MUBs in low dimensions 

In what follows we will assume that a complete system of MUHs 
Hx, ■ ■ ■ Hd is given. In fact, much of the discussion below remains valid 
for non-complete systems after appropriate modifications, but it will 
be technically easier to restrict ourselves to the complete case. The 
general aim is to establish structural properties of Hx, . . . Hd which give 
restrictions on what a complete system may look like. If some of these 
properties were to contradict each other in a non-primepower dimension 
d, then we could conclude that a complete system of dimension d does 
not exist. This is one of the main tasks for future research, mainly 
for d = 6. We will give some non-existence results in this direction in 
Section [3j 

Consider each appearing complex Hadamard matrix Hj as a d-element 
set in T d (the elements are the columns Ci, . . . c<j of the matrix; the de- 
pendence on j is suppressed for simplicity), and introduce its Fourier 
transform 

d 

(1) g^) := 4-( 7 ) = c l for each 7 <= Z d . 

k=X 

Notice that the orthogonality of the rows of Hj implies that if p G Z d 
is any permutation of the vector (1, —1, 0, 0, . . . , 0) then 

(2) 9j(p)=0- 

Also, note that conjugation is the same as taking reciprocal for uni- 
modular numbers, i.e. ^(7) = Ylt=i c fc 7 > an< ^ therefore the square of 
the modulus of 5^(7) can be written as 

d 

(3) 0,(7) := M7)| 2 = £>*/ c l) 7 for each 7 G Z d . 

k,l=X 
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Also, introduce the notation 

d 

(4) G(rf) := G j( r y) for each 7 G z<i - 

In similar fashion, introduce the Fourier transform of the whole system 
as 

d 

(5) /(7) := Xl^'W foreach 7 GZ d , and 

i=i 

(6) F( 7 ) := |/( 7 )| 2 = E( u / y ) 7 for each ^ e ^ 

u,v 

where the summation goes for all pair of columns u, v in the matrices 
H\, . . . , H d . 

The main advantage of taking Fourier transforms is that any polyno- 
mial relation (such as orthogonality or unbiasedness) among the entries 
of the matrices Hj will be turned into a linear relation on the Fourier 
side. We will collect here linear equalities and inequalities concerning 
the functions .F( 7 ) and G("f). 

Let 7r r = (0, 0, ... 0, 1, 0, ... 0) G Z d denote the vector with the rth 
coordinate equal to 1. Then for each j — 1, . . . d we have 

d d / d \ d f d \ 

$>,(7+7r r ) = 5>*/ c <) 7+7rr = E( c * /c ') 7 E^/^)" 

r=l r=l \k,l=l J k,l=l \r=l J 

and observe that the last sum is zero by orthogonality if k 7^ I, while 
it is d if k — I. This means that for each j — 1, . . . d, 

d 

( 7 ) Yl G ^ + K r) = ^ fOT 6aCh 7 G Z ^ 
r=l 

which then implies 

d 

(8) G( 7 + 7r r ) = rf 3 for each 7 G Z d . 

r=l 

In a similar fashion we can turn the unbiasedness relations also to 
linear constraints on the Fourier side. Let u/v = (z±, z-i . . . , Zd) G T d 
be the coordinate-wise quotient of any two columns from two different 
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matrices from Hi, . . .Ha- Then u and v are unbiased, which means 
that 

(9) = | J^Zr| 2 ~d = ^Zr/Zt. 

r r^t 

Using this we can write 

(10) £(F - G)(7 + vr r - 7T t ) = ^(u/v)^ ( ^(u/v)— ) = 0, 

r^t u,v \ r^t / 

where the summation on u, v goes for all pairs of columns from different 
matrices, and the last equality is satisfied because each inner sum is 
zero by (EJ). Also, by (ED we have dG{j) + J2 r & G< CI + - *t) = d 4 , 
and we can use this to rewrite ffTOl as 

(11) dG( 1 ) + J2 F ^ + 7Y r-n t ) = d\ 

which is somewhat more convenient than ( |T0l ). 

We also have some further trivial constraints on F and G. Namely, 

(12) F(0) = d\ G(0) = d 3 , and 

(13) < F( 7 ) < d\ < G(i) < d 3 , for each 7 G Z d . 

Also, by the inequality of the arithmetic and quadratic means we have 

(14) F( 7 ) < dG(j), for each 7 e Z d . 

The point is that the linear constraints flHj). fTTTj) . (fT2j). f[T3]) . (1T4D put 
severe restrictions on the functions F and G. In fact, it turns out that 
all the structural results on complete systems of MUHs in dimensions 
2, 3, 4, 5 follow from these constraints. These structural results are not 
new (cf. [5]) but nevertheless we list here the two most important ones 
as an illustration of the power of this Fourier approach. The first one is 
a celebrated theorem of Haagerup [IDj which gives a full classification 
of complex Hadamard matrices of order 5. In the original paper \X0\ 
the author combines several clever ideas with lengthy calculations to 
derive the result, whereas it follows almost for free from the formalism 
above. 

Proposition 2.1. Any complex Hadamard matrix of order 5 is equiv- 
alent to the Fourier matrix F5, given by F^{j, k) = u^" 1 '^" 1 ' , (j, k = 
!,...,§), where u = e 2ln ^ 5 . 
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Proof. Let Hi be a complex Hadamard matrix of order 5. Then the 
function G\{pf) = | -^1 (t) 1 2 satisfies equation (J7J) for all 7 £ Z 5 , and we 
have Gi(0) = 25 and < Gi(j) < 25 for all 7 £ Z 5 . Regarding each 
G\ (7) as a nonnegative variable (as 7 ranges through a sufficiently large 
cube around the origin in Z 5 ), a short linear programming code testifies 
that under these conditions G\(p) = 25 for all such p £ Z 5 which is 
a permutation of (5, —5, 0, 0, 0). Also, we may assume without loss of 
generality that Hi is normalized (i.e. its first row and column are made 
up of Is), and then the information G\(p) = 25 implies that all other 
entries of Hi are 5th roots of unity. It is then trivial to check that there 
is only one way (up to equivalence) to build up a complex Hadamard 
matrix from 5th roots of unity, namely the matrix F 5 . □ 

We remark here all the linear programming mentioned in this paper 
uses rational coefficients, so no numerical errors are encountered, and 
each result is certifiable. Let us also remark that Proposition 12.11 is 
the only non-trivial result concerning MUHs and MUBs in dimensions 
d < 5. The classification of complex Hamamard matrices and MUBs is 
more or less trivial for d — 2, 3, 4 due to the geometry of complex unit 
vectors. We give here the essence of this classification (for full details 
see [5]). 

Proposition 2.2. In any normalized complete system of MUHs in di- 
mension d = 3,4, 5 all entries of the matrices are dth roots of unity. 
For d = 2 all entries are 4th roots of unity. 

Proof. Let d = 3,4,5. Assume Hi, . . . Hd is a normalized complete 
system of MUHs. Then the functions F and G must satisfy the linear 
constraints ©, (HH), (IH, AH; flU- Regarding each F( 7 ) and G(j) 
as a nonnegative variable (as 7 ranges through a sufficiently large cube 
around the origin in Z d ), a short linear programming code testifies 
that under these conditions F(p) = d 4 for all such p £ Z d which is a 
permutation of (d, —d, 0, . . . , 0). This means that all entries in all of 
the matrices must be dth roots of unity. The proof is analogous for 
d = 2 except that in this case we can only conclude F(4, —4) = 16, so 
that the matrices contain 4th roots of unity. □ 

Let us make a remark here about d = 4. In this case it is not true 
that all normalized Hadamard matrices must be composed of 4th roots 
of unity. However, it is true that a complete system of MUHs must 
be composed of such. This phenomenon shows up very clearly in our 
linear programming codes. Writing the constraints ([7]) on (^1(7), and 
L?i(0) = 16, and < (^1(7) < 16 does not enable us to conclude that 
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Gi(p) = 16 with p being a permutation of (4, —4, 0, 0. However, writing 
all the constraints (fSj). (ITTj). (fT3|). (113]), (TT4|) on the functions F and G 
we can indeed conclude that F(p) = 4G(p) = 256. 

We end this section with a few remarks concerning d = 6. If we 
could similarly conclude that 

(15) F(p) = 6 4 for all p being a permutation of (6, —6, 0, 0, 0, 0) 

then it would mean that a complete system of normalized MUHs in di- 
mension 6 can only be composed of 6th roots of unity. Such a structural 
information would be wonderful, as it is easy to check by computer that 
no such complete system of MUHs exists. Therefore, we could conclude 
that a complete system of MUHs does not exist at all. Unfortunately, 
the constraints ©, (HU), (H3), (USD, (HID do not seem to imply flTS]). At 
least, we have run a linear programming code with 7 ranging through as 
large a cube as possible (due to computational limitations), and could 
not conclude ffToT) . Nevertheless, our main strategy for future research 
in dimension 6 must be as follows: using the linear constraints on F 
and G try to establish some structural information on the vectors ap- 
pearing in a hypothetical complete system of MUHs, and then show by 
other means (e.g. a brute force computer search) that such constraints 
cannot be satisfied. We formulate here one conjecture which could be 
crucial in proving the non-existence of a complete system of MUHs in 
dimension 6. 

Conjecture 2.3. Let Hi be any complex Hadamard matrix of order 6, 
not equivalent to the isolated matrix Sq (cf. [T7j for the matrix Sq). Let 
p be any permutation of the vector (1, 1, 1, —1, —1, —1). Then gi(p) = 
for the function g± defined in ([ID . 

This conjecture is supported heavily by numerical data. We have 
tried hundreds of matrices randomly from each known family of com- 
plex Hadamard matrices of order 6 (a complete classification is unfor- 
tunately not available). Currently we cannot prove this conjecture, but 
in Section 3 we will show an example of how it could be used in the 
proof of non-existence results (cf. Remark [3.41) . 

3. Non-existence results 

We now turn to non-existence results, namely that complete systems 
of MUHs with certain properties do not exist. The first of these, which 
we regard as the main result of the paper, is that any pair of real 
unbiased Hadamard matrices cannot be part of a complete system of 
MUHs. In fact, we prove the following stronger statement. 
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Theorem 3.1. Let Hi, . . . Hd be a complete system of MUHs such that 
Hi is a real Hadamard matrix. Then any column vector v = (vi, . . . , Vd) 
of the other matrices H 2 , . . . Hd satisfies that Ylt=i v 1 = 0- 

Proof. Let 7^ p = (r 1 ,...,r d ) G Z d be such that J2k=i r k = 
and J2t=i \ r k\ — 4. There are five types of these vectors (up to 
permutation): (1, -1, 0, . . . , 0), (2, -2, 0, . . . , 0), (2, -1, -1, 0, . . . , 0), 
(-2, 1, 1, 0, . . . , 0), and (1, 1, -1, -1, 0, ... 0). Then, Theorem 8 in 
(or Corollary 2.4 in [16J) shows that the function / defined in ([5]) sat- 
isfies 

(16) f(p) = 

for all these vectors p. 

Let Ci, c 2 , . . . , c^2 denote the column vectors appearing in the system 
H u ...H d . For each 7 G Z d let 

(17) v( 7 ) = (cL...c] 2 )GT d2 

for k = l,...d. Consider the vectors 7 fc = (0, ... 0, 2, 0, ... 0) G 7L d 
with the 2 appearing in position k. Finally, consider the vector w = 
^fe=i v (7fc), and let us evaluate ||w|| 2 . On the one hand, the vectors 
v(7fc) are all orthogonal to each other by ( fl~6l) . and they all have length 
ll v (7fc)l| 2 = anc l hence ||w|| 2 = d 3 . On the other hand we know 
the first d coordinates of w. Each v( 7fc ) has first d coordinates equal 
to 1, because Hi is a real Hadamard matrix. Therefore the first d 
coordinates of w are all equal to d. Therefore, ||w|| 2 > d 3 on account 
of the first d coordinates. Hence, all other coordinates of w must be 
zero, which is exactly the statement of the theorem. □ 

Theorem 13.11 implies immediately the following corollary. 

Corollary 3.2. Let Hi,...Hd be a complete system of MUHs such 
that Hi is a real Hadamard matrix. Then there is no further purely 
real column in any of the matrices H 2 , . . . , H^- In particular, it is 
impossible to have two real Hadamard matrices in a complete set of 
MUHs. 

This statement is sharp in the sense that for d = 2, 4 the complete 
systems of MUHs are known to contain one real Hadamard matrix. 
Also, in several dimensions d = An 2 pairs (and even larger systems) of 
real unbiased Hadamard matrices are known to exist [U [11], so that 
the corollary above is meaningful and non-trivial. 

Our next result is a new proof of the fact in dimension 6 the Fourier 
matrix F e cannot be part of a complete system of MUHs. This result 
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is well-known, but the only proof we are aware of uses some computer 
algebra, while we present an easy conceptual proof here. 

Proposition 3.3. There exists no complete system of MUHs in di- 
mension 6 which contains the Fourier matrix F$ . 

Proof. Assume by contradiction that such a system Hi, . . . H 6 exists, 
and assume Hi = F 6 . Consider the vectors 71 = (1,1,1,0,0,1), 
72 = (0,0,1,1,1,1), 73 = (1,1,0,1,1,0), 7 4 = (0,1,0,1,0,2), 75 = 
(1, 0, 0, 0, 1, 2), and 76 = (0, 1, 0, 0, 2, 1), and consider the corresponding 
vectors v(7 fc ) defined in fTl7|) . and let w = Ylt=i v (7fc)- All the vectors 
v(t&) are orthogonal to each other by ( 116ft . therefore ||w|| 2 = 216. On 
the other hand, we know the first 6 coordinates of w. It is easy to 
calculate that each of these coordinates has modulus 6, and therefore 
||w|| 2 > 216 on account of the first 6 coordinates. This implies that 
all the other coordinates of w must be zero. This yields a polynomial 
identity for the coordinates of any column vector appearing in the ma- 
trices if 2, • • • , Hq. Instead of using this identity directly, however, we 
observe that the same argument applies to the vectors 71, . . .75 and 
7e = (2, 0, 0, 1, 0, 1), and w' = v(t£) + Ylt=i v (lk)- By considering the 
difference w — w' we conclude that v(7@) and v(t£) must coincide in 
the last 30 coordinates. That is, if (z±, . . . , Zq) is any column vector in 
the matrices H 2 , . . . H e then z 2 z\z% = zfz^z^, and hence z 2 z\ = z\z±. 
Furthermore, one can permute the coordinates of 7^ in a cyclic manner, 
and the argument remains unchanged, yielding this time z§z\ = z\zi. 
Dividing these two equations finally gives z 5 /z 2 = z\jz± for each of the 
last 30 vectors in our complete system of MUHs. This means, by defi- 
nition, that the last 30 coordinates of the vectors v(0, —1, 0, 0, 1, 0) and 
v(l, 0, 0, — 1, 0, 0) coincide. But this is a contradiction, because these 
vectors should be orthogonal to each other by ([161) . □ 

Remark 3.4. Finally, we discuss informally a non-existence result in 
which we use Conjecture 12.31 in the proof. Nevertheless, the result 
itself is not "conditional" because it was proved earlier in [13] by a 
massive computer search after a discretization scheme. The argument 
we present here is much more elegant though, and shows a possible way 
forward in proving the non-existence of complete systems of MUHs in 
dimension 6. 

We claim that there exists no complete system Hi, ... , H$ of MUHs 
in dimension 6 which contains any of the matrices F§{a, b) of the Fourier 
family (cf. [IT] for the Fourier family Fe(a, b)). We sketch the proof 
here, on the condition that Conjecture 12.31 is valid. First of all, the 
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statement is equivalent for the transposed family Fjf(a,b), as it is 
well-known that a matrix H can be part of a complete system of 
MUHs if and only if H T can. The significance of this is that each 
member of the transposed family Fq(cl, b) contains the three vectors 
Cx = (l,l,l,l,l,l),c 2 = {1,uj,uj 2 ,1,uj, co 2 ) and c 3 = (l,u 2 ,u,l, co 2 , co), 
where co = e 2l7T ^ 3 . Also, it is well-known that a complex Hadamard ma- 
trix equivalent to Sq cannot be part of a complete system of MUHs, so 
that we can assume without loss of generality that Hi, . . . Hq are not 
equivalent to Sq. The significance of this fact is that now Conjecture 
12.31 (if true) can be invoked. One can make a clever selection of vectors 
7ij ■ ■ -7i2 £ 7$ such that the same argument as in Proposition 13.31 can 
be used. Namely, all the vectors v(7fc) are orthogonal to each other by 

either (USD or by Conjecture [231 so that ||w|| 2 = | Y%=\ v (^)l 2 = 432 - 
(This is where we use Conjecture 12.31 otherwise the appropriate se- 
lection of 71, . . . 712 would not be possible.) On the other hand, three 
coordinates of w corresponding to the columns Ci, C2, C3 are known ex- 
actly, and the modulus of them happens to be 12, again by the choice 
of the vectors 7^. This leads us to conclude that all the other 33 co- 
ordinates of w must be zero, which yields a polynomial identity for 
the coordinates of each the 33 unknown columns in Hi, . . . , H e . Fi- 
nally, one can make several such selection of 71, . . . 712, each yielding a 
polynomial constraint on the unknown columns, and these constraints 
contradict each other just as in the proof of Proposition 13.31 □ 

We believe that the proof of the non-existence of complete systems 
of MUHs in dimension 6 will hinge on Conjecture 12.31 The reason 
is that it introduces yet another non-trivial linear constraint on the 
function G, and these constraints will ultimately lead to a contradiction 
(maybe indirectly, as in Proposition 13. 3p . Therefore, we would be very 
interested to see a proof of Conjecture 12.31 
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